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We study finite-temperature phase transition and equation of state for two-flavor QCD at TVt = 4 using an 
RG-improved gauge action and a meanfield-improved clover quark action. The pressure is computed using the 
integral method. The 0(4) scaling of chiral order parameter is also examined. 



1. Introduction 



The transition temperature and the equation of 
state (EOS) of QCD at finite temperature belong 
to the most basic information for understanding 
the early Universe and heavy ion collisions. Full 
QCD studies of these quantities have been made 
mainly with the Kogut-Susskind quarks, particu- 
larly for the EOS In this paper we present the 
first result of the EOS from Wilson-type quarks. 

We study two-flavor QCD on A^t = 4 lattices. 
In order to suppress lattice artifacts, which are 
known to be severe for the combination of the pla- 
quette gauge and Wilson quark actions, we adopt 
a renormalization-group (RG) improved gauge 
action |^ combined with a meanfield-improved 
clover quark action. See Ref. Q for details of 
our action. 

In Fig. |l] we show the phase diagram with our 
action at Nt — 4. The line of finite-temperature 
transition is determined by the Polyakov loop and 
its susceptibility. The parity-broken phase Q is 
not yet identified. Dashed lines are used in a test 
discussed in Sec. 0. 



0.05 




O— OK, 

• data 1 6^x4 

path |3=2.10 

path |3=2.20 

path K=0.13 



Figure 1. Phase diagram for Nt = 4. The solid 
line represents the critical line KdP) of vanish- 
ing pion mass at T = 0. The shaded line Kt{/3) 
is the location of finite-temperature transition. 
Dots represent the simulation points on 16'^ x 4 
and 16^ lattices carried out so far in the present 
work. 
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Figure 2. Pressure on an 8^ x 4 lattice calculated Figure 3. Pressure on a 16^ x 4 lattice as a 

from three paths in Fig. ^ function of tops/tovIt=o- 



2. Equation of state 

We compute the pressure p by the integral 
method which is based on the formula, valid 
for large homogeneous systems, that 

The integration path in the parameter space 
(/3, K) should start from a point in the low 
temperature phase where the integrand approxi- 
mately vanishes. We evaluate the quark contribu- 
tions to the derivatives ^ and ^ by the method 
of noisy source using U(l) noise vectors. 

The value for the pressure computed in (|l|) 
should be independent of the choice of the in- 
tegration path. To check this point, we make a 
series of test runs on S'^ x 4 and S** lattices along 
three paths shown in Fig. ^ generating 500 HMC 
trajectories at each point. The results for p/T'*' 
obtained from these paths are plotted in Fig. 
We find that p/T^ at {/3,K) = (2.1,0.13) and 
(2.2,0.13) in the two figures are in good agree- 
ment, confirming the path independence of the 
integral. 

Encouraged by this result, we perform produc- 
tion runs on 16^ x 4 and 16'* lattices. At each 
dots plotted in Fig. |l|, we generate 500-2000 tra- 
jectories on the 16^ x 4 lattice and 200-300 tra- 
jectories on the 16^ lattice. Measurement of the 



derivatives is made at every trajectory. Hadron 
propagators are calculated at every fifth trajec- 
tory to compute pseudo scalar (mps) and vector 
(mv) meson masses. 

As is seen from Fig. ||, paths along the in- 
direction give much smaller errors in p/T" than 
those from paths in the /3-direction. Therefore, 
we carry out the integral in the i^-direction. We 
then obtain pressure plotted in Fig. ^ as a func- 
tion of the mass ratio mps /my at zero tempera- 
ture. Interpolating these data, we find p/T^ for 
each value of mps/my. 

Figure ^ shows the pressure as a function 
of T/Tpc at fixed mps/rny. Here Tpc is the 
pseudo-critical temperature at the same value of 
"^ps/"^v■ The temperature scale is set by my 
through T/Tpc = rny (Ppc) / my (P) with l3pc the 
pseudo-critical coupling. 

We find that the pressure for fixed T/Tpc de- 
pends only weakly on the quark mass even for rel- 
atively heavy quark in the range mps/my = 0.7- 
0.8. For heavier quark masses, the pressure de- 
creases toward the pure gauge value (dashed line) 
P as expected. 

We also note that the magnitude of pressure is 
much larger than that for the pure gauge system 
for Nt — 4, and that it overshoots the Stefan- 
Boltzman value in the continuum at high temper- 
atures. These features are probably the result of 
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Figure 4. Pressure on a 16"^ x 4 lattice as a 
function of T /Tpc- The dashed curve shows pres- 
sure for pure gauge theory with the RG-improved 
action on a 16^ x 4 lattice [p|. 



large discretization errors from the clover quark 
action here . Indeed the large Stefan-Boltzman 
value on an A'^f = 4 lattice shown at the top-right 
in Fig. ^ is dominated by the quark contribution. 

3. 0(4) Scaling 

The chiral phase transition of two-flavor QCD 
is expected to belong to the universality class of 
0(4) spin system in three dimensions. In par- 
ticular, identifying the magnetization, external 
magnetic field, and reduced temperature of the 
spin model with M = (^5'), h = 2mqa, and 
t = j3 — /3ct, where Pet is the chiral transition 
point, we expect a scaling relation 



(2) 



to hold with the 0(4) scaling function f(x) ^ and 
the 0(4) critical exponents (3 and (5 0. A pre- 
vious study using the RG-improved gauge action 
and unimproved Wilson quark action ||7|] found 
this relation to be well satisfied for the quark mass 
and the chiral condensate defined by axial Ward 
identities ||ic| ] . 

Figure |5| shows the result of a similar analy- 
sis from the present work. Data for (^"I'}sub ~ 
2mqa(2X)^ ^^(7r(a;)7r(0)) are fitted to the scal- 
ing relation, adjusting Pet and the scales for t 
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Figure 5. 0(4) scaling relation. 



and h. The scaling ansatz works well, yield- 
ing Pet = 1.47(7) for the best fit of data for 
2m,a < 0.9 and P < 1.95 with xVdf : 
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